Non-Fermi liquid behavior with and without quantum criticality in Cei_ x Yb z CoIn5 
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One of the greatest challenges to Landau's 
Fermi liquid theory - the standard theory of 
metals - is presented by complex materials with 
strong electronic correlations. In these materials, 
non-Fermi liquid transport and thermodynamic 
properties are often explained by the presence 
of a continuous quantum phase transition which 
happens at a quantum critical point (QCP). A 
QCP can be revealed by applying pressure, mag- 
netic field, or changing the chemical composi- 
tion. In the heavy- fermion compound CeCoIns, 
the QCP is assumed to play a decisive role in 
defining the microscopic structure of both normal 
and superconducting states. However, the ques- 
tion of whether QCP must be present in the ma- 
terial's phase diagram to induce non-Fermi liquid 
behavior and trigger superconductivity remains 
open. Here we show that the full suppression of 
the field-induced QCP in CeCoIns by doping with 
Yb has surprisingly little impact on both uncon- 
ventional superconductivity and non-Fermi liquid 
behavior. This implies that the non-Fermi liquid 
metallic behavior could be a new state of matter 
in its own right rather then a consequence of the 
underlying quantum phase transition. 

The heavy-fermion material CeCoIns is a prototypi- 
cal system in which strong interactions between conduc- 
tion and predominantly localized /-electrons give rise to a 
number of remarkable physical phenomena [TJ|2]- Uncon- 
ventional superconductivity emerges in CeCoIn 5 out of a 
metallic state with non-Fermi-liquid (NFL) properties: 
linear temperature dependence of resistivity below 20 K, 
logarithmic temperature dependence of the Sommerfeld 
coefficient, and divergence of low temperature magnetic 
susceptibility [3HS]- These anomalies disappear beyond 
a critical value of the magnetic field and the system re- 
covers its Fermi liquid properties. The crossover from 
non- Fermi liquid to Fermi liquid behavior is thought to be 
governed by a quantum critical point (QCP), which sepa- 
rates paramagnetic and antiferromagnetic (AFM) phases 
and is located in the superconducting phase [7J [5] ■ Neu- 
tron scattering studies [S] and more recent measurements 
of the vortex-core dissipation through current-voltage 
characteristics [10 provide direct evidence for an anti- 
ferromagnetic QCP in CeCoIns that could be accessed 
by tuning the system via magnetic field or pressure. 



Nevertheless, a growing number of /-electron systems 
do not conform with this QCP scenario; for example, the 
NFL behavior in some systems occurs in the absence of 
an obvious QCP [11] [12]. An intriguing candidate is Yb- 
doped CeCoIns that exhibits an unconventional T — x 
phase diagram without an apparent QCP, while the on- 
set of coherence in the Kondo lattice and the supercon- 
ducting transition temperature T c are only weakly de- 
pendent on Yb concentration and prevail for doping up 
to x = 0.65. [13]. Yet, the presence of a QCP in the 
parent CeCoIns compound and the logarithmic temper- 
ature dependence of normal state Sommerfeld coefficient 
in lightly doped Cei-^Yb^CoIns crystals [H] show that 
this system is in the vicinity to a QCP. Therefore, it is 
important to elucidate the nature of the NFL behavior 
and unconventional superconductivity in such a system, 
to search for possible QCPs, and to determine the degree 
to which quantum criticality and superconductivity are 
coupled to each other. 

Here we focus on revealing a field-induced quantum 
critical point through normal-state magnetoresistivity 
measurements and finding its evolution with Yb-doping. 
Experimental details are given in the "Materials and 
Methods" Section. We first present the results of our 
study of the magnetic field H and temperature depen- 
dence of the transverse (H _L ab) in-plane magnetore- 
sistivity Api/p a = [p^{H) - p o (0)]/p„(0), for H < 14 
T and 3 < T < 70 K. Figure 1(a) and its inset dis- 
play the field dependence of Ap^ / p a measured at differ- 
ent temperatures for the x = and x = 0.4 samples, 
respectively. We note that the data for these samples 
fall into two groups: (i) non-monotonic field dependence 
of magnetoresistivity (MR) [Fig. 1(a)] with quadratic 
MR at high fields, typical for x < 0.20 and (ii) negative 
and quadratic MR over the whole measured field range 
[inset to Fig. 1(a)], typical behavior for the high Yb 
(0.25 < x < 0.65). We, therefore, conclude that the MR 
for low Yb doping has two main contributions: one nega- 
tive and quadratic in H, which we denote Ap° rb / p a [black 
symbols of Fig. 1(b) at high fields] and the other is posi- 
tive denoted Ap s P m / p a [green symbols of Fig. 1(b)], with 
the latter one obtained by subtracting for all field val- 
ues the negative quadratic MR from the measured MR. 
This latter contribution to MR is isotropic since the green 
data follow quite well the longitudinal magnetoresistivity 
Apl/p a [red data of Fig. 1(b)]. Also, Ap spm / p a is linear 
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FIG. 1: (Color online) (a) Magnetic field scans of trans- 
verse magnetoresistivity Ap^ / p a measured at different tem- 
peratures on Cei-^Yb^CoIns for x — (main panel) and 
x = 0.4 (inset), (b) In-plane transverse Ap^/p a (black) and 
longitudinal Ap\ / p a (red) magnetoresistivities (MR) of x = 
and x — 0.4 samples vs H 2 measured at 5 K. The green data 
are denoted Ap apm lp a and represent the component of MR 
obtained by subtracting the high field quadratic MR from the 
measured transverse MR. (c) Apa pl "/p a vs H (for low H val- 
ues) for the x = sample measured at different temperatures, 
(d) H dependence of the isotropic component Ap spm I p a of 
transverse magnetoresistivity of CeCoIn 5 measured at differ- 
ent temperatures and up to 14 T. 



in H at low fields [Fig. 1(c)] and saturates at high fields 
[Fig. 1(d)]. 

Positive MR of heavy-fcrmion materials at low fields 
marks the departure from the single-ion Kondo behav- 
ior and is determined by the formation of the coherent 
Kondo lattice state in systems in or close to their Fermi- 
liquid ground state [T5H20] . The maximum in the MR 
of a Kondo-lattice Fermi liquid at a certain field value 
is a result of the competition between a T-independent 
residual resistivity contribution that increases with in- 
creasing H , and a T-dependent term that decreases with 
increasing H [20j[2T]. Thus, in conventional Kondo lat- 
tice systems, the peak in the field-dependent MR moves 
toward lower H with increasing T since a lower field is 
required to break the Kondo lattice coherence. 

To determine the nature of the positive magnetore- 
sistance in Cei-zYb^CoIns for x < 0.20, we extract 
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FIG. 2: Temperature dependence of the characteristic field 
H max corresponding to the maximum isotropic component of 
magnetoresistance MR^™. The solid line is a guide to the 
eye, while the dotted line is a linear fit of the low-T linear 
behavior of H max for the x = sample. 



the temperature dependence of the field H max at which 
Apa Pm /p a saturates at its maximum value Mi?^™ [see 
Fig. 1(d)]. The data of Fig. 2 show that the temper- 
ature dependence of H max of the x = sample is non- 
monotonic, with a maximum around 20 K, and a linear 
behavior at low-T values. This non-monotonic T behav- 
ior of H rnax is representative for all the samples with 
x < 0.20. The positive MR measured at T > 20 K [for 
which H max (T) decreases with increasing T] could re- 
flect the presence of the coherent Kondo lattice state at 
low field values as discussed in the previous paragraph. 
In contrast, the behavior below 20 K is opposite with 
the one discussed above for conventional Kondo lattice 
systems. The increase of H max with increasing T at 
these lower temperatures had previously been observed 
in CeCoIns and it has been attributed to field quenching 
of the AFM spin fluctuations responsible for the NFL be- 
havior [22 . Therefore, we conclude that the positive MR 
measured at T < 20 K in Cei-^Yb^CoI^ with x < 0.20 
reflects the dominant role played by the AFM quantum 
spin fluctuations. 

An important next goal is to identify the quantum crit- 
ical field (Hqcp) associated with these quantum fluctu- 
ations. One option is to extrapolate the low tempera- 
ture linear H max (T) behavior to T = K and identify 
this H max (0) with Hqcp- However, since there is a cer- 
tain error associated with the determination of H max , we 
adopted a more accurate procedure to unambiguously de- 
termine Hqcp for different Yb doping. This procedure is 
discussed in detail in the "Materials and Methods" Sec- 
tion below. We show in Fig. 3(a) the values of Hqcp as 
a function of Yb concentration. As expected, the value 
of Hqcp of 4.1 T for CeCoIn 5 coincides with the value of 
Hqcp determined previously from both resistivity mea- 
surements done in the normal state 1 23 and I — V char- 
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FIG. 3: (Color online) (a): Evolution of the field induced 
quantum critical point Hqcp, coherence temperature T co h, 
superconducting critical temperature T c (from Ref. [13]) of 
Cei- x R x Coln^ (R =Yb, La) as a function of rear-earth con- 
centration x. (b): Plot of T co h vs. T c for Yb- and La-doped 
(from Ref. [3]) samples. 



acteristics measured in the mixed state [TU]. Therefore, 
the measurement of Ap^ / p a along with the analysis used 
here represent an excellent experimental technique to de- 
termine the field-induced QCP in the NFL regime. The 
Hqcp = points in Fig. 3(a) for x > 0.20 correspond to 
the case when there is no positive MR, i.e., the maximum 
in MR shifts to zero field [see inset to Fig. 1(a)]. 

Figure 3(a) also shows the suppression of T c (from [13] ) 
and T co h with doping for Ce 1 _ a ,Yb 2 .CoIn 5 . These doping- 
dependent Hqc p and temperature phase diagrams show 
that while superconductivity (SC) is robust and survives 
over the whole Yb doping range, the field-induced QCP 
is strongly suppressed with Yb doping and disappears 
for x > 0.20. This implies that SC and quantum criti- 
cality are likely to be decoupled in this system, i.e., un- 
conventional superconductivity is not triggered by spin 
fluctuations. 

The experimental technique used to determine Hqcp 
also permits the determination of the gyromagnetic fac- 
tor g (see "Materials and Methods" Section below). There 
is a significant change in the value of the g factor at QCP, 
from 2.2 just below Hqcp to 1.3 just above Hqcp- This 
means that for H < Hqcp the conduction electrons be- 
come weakly coupled to the local spins, so that the heavy- 
fermions would have a lower effective mass. However, for 
H > Hqcp the antiferromagnetic fluctuations between 
the local moments are suppressed, so that the heavy- 
fermions recover their effective mass and this is reflected 
in the reduction in the value of the g- factor |25] . 

Figure 3(b) shows that T co h (from this work) scales 
with T c (from Ref. [13]) over the whole Yb doping stud- 
ied here (x < 0.65) (black squares) and La doping (red 
circles, data from Ref. [3])- We note that T coh vs T c fol- 
lows a linear dependence for the whole Yb doping with 
x < 0.65 and up to about 10% La doping, with the two 



slopes having about the same value; this figure also shows 
the result of the linear fit for the Yb doping. Thus, it is 
intriguing that while the scaling of these two tempera- 
tures is in conformity to the other rare-earth substitu- 
tions on the Ce site [23], the nature of the scaling sets 
the Yb substitution aside from the other rare-earth sub- 
stitutions since, in the latter ones a non linear and quite 
fast suppression of both temperatures is observed, with 
a suppression of T c to zero at around x = 17.6% of rare- 
earth substitution (see Fig. 3(a), data from [31 [T3]). 

A very interesting and puzzling behavior of the 
Cei_ x Yb x CoIn 5 system is that, even though QCP disap- 
pears for x > 0.20, the system continues to display NFL 
behavior, as evidenced by the sublinear T-dependence of 
its resistivity [13] . This means that this non- Fermi liquid 
behavior at higher Yb doping could be a new state of 
matter in its own right rather then a consequence of the 
underlying quantum phase transition. We further inves- 
tigated the origin of this NFL behavior by studying in 
more detail the T dependence of the resistivity measured 
in different magnetic fields. The resistivity of all the Yb- 
doped single crystals studied follow remarkably well the 
expression Pa(H) = po + AT + By/T for temperatures 
up to about 15 K [see inset to Fig. 4(a) for the fitting 
results] and in fields up to 14 T; here po, A, and B are 
doping- and field-dependent fitting parameters. We show 
in Figs. 4(a) and 4(b) the doping and magnetic field de- 
pendences, respectively, of the ratio of the T-linear con- 
tribution to the total T-dependent contribution of the 
resistivity. These results show that the linear in T term 
in resistivity is present for low Yb doping (x < 0.20) in 
the quantum critical regime (T < 20 K). The percentage 
of the linear in T term decreases with increasing x and 
H and disappears for x > 0.20 and at fields at which the 
MR is only negative, where only the y/T dependence is 
present. The bottom inset to Fig. 4(a) shows the av- 
erage Yb valence as a function of Yb doping [2B] . It Is 
noteworthy that the average Yb valence decreases with 
increasing Yb and saturates to a value of about 2.3 for 
x > 0.20 [SHUT]. Th is result along with the data of Figs. 
4(a) and 4(b) and the fact that both the superconduct- 
ing transition and coherence temperatures remain weakly 
dependent on doping indicate that Yb atoms form a co- 
operative mixed-valence state that significantly reduces 
the pair-breaking effects, which could also play an im- 
portant role in the origin of the NFL behavior at these 
higher Yb concentrations (x > 0.20). This idea is sup- 
ported by the observed linear dependence of T co ^ vs T c 
[Fig. 3(b)]. The "Materials and Methods" Section gives 
more discussion and interpretation of this scaling. 

All of the above observations show that the field- 
induced QCP plays a supporting rather then competing 
role in the emergence of unconventional superconductiv- 
ity. Moreover, we conclude that a highly unusual mi- 
croscopic mechanism of unconventional superconductiv- 
ity in Cei-^Yb-rCoIns may be at play: on one hand, the 
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FIG. 4: (color online) (a) Doping dependence of linear in 
temperature T contribution devided by the total contribu- 
tion to the resistivity obtained from fits (shown in the top 
inset) of the resistivity data by p^ (H) — po + AT + By/T, 
along with the evolution of the field-induced QCP Hqcp(x) 
(main panel). Bottom inset: Variation of average Yb valence 
with doping (data from Ref. 26 ). (b) Field dependence of 
the linear in T contribution over the total contribution to 
the resistivity (left axis) and magnetoresistivity (right axis) 
measured at 3 K for pure CeCohis. 



normal state is fully reconstructed by the ytterbium sub- 
stitution, while on the other hand, the superconducting 
critical temperature is reduced only by a factor of two 
at x ~ 65%, strongly suggesting that the superconduct- 
ing pairing is spatially inhomogeneous and involves local 
Ce /-moments. Finally, we propose a novel technique to 
probe the interplay between quantum criticality and su- 
perconductivity, which can be used to analyze a variety 
of strongly correlated electronic materials. 



the field-induced quantum critical point from magnetore- 
sistivity data. Finally, we offer an interpretation of the 
scaling of the Kondo lattice coherence temperature and 
superconducting critical temperature. 



Experimental methods 

Single crystals of Cei-^Yb^CoIns (0 < x < 0.70) were 
grown using an indium self - flux method |28j . The crys- 
tal structure and composition were determined from X- 
ray powder diffraction (XRD) and energy dispersive X- 
ray (EDX) techniques. The single crystals have a typical 
size of 2.1 x 1.0 x 0.16 mm 3 , with the c axis along the 
shortest dimension of the crystals. They were etched 
in concentrated HC1 for several hours to remove the in- 
dium left on the surface during the growth process and 
were then rinsed thoroughly in ethanol. Four leads were 
attached to the single crystal, with I \\ a. High qual- 
ity crystals were chosen to perform in-plane transverse 
(Apa/p a ) and longitudinal {Ap\/ p a ) magneto-resistivity 
(MR) measurements, with H J_ ab and H \\ ab, respec- 
tively, as a function of temperature (T) and applied mag- 
netic field (H). In both situations, however, the field 
is perpendicular to the current, H J_ /, to ensure that 
Lorentz force remains the same. 



Determination of field-induced quantum critical 
point 

We used to following procedure to determine the quan- 
tum critical field Hqcp for different Yb doping. We 
define the characteristic fields H max where Ap* pm / ' p a 
is 100%, 95%, 80%, 75%, etc. of the saturation value 
MRfP™. We show the T dependence of these character- 
istic fields in Fig. 5(a) by the black, red, green, and blue 
solid circles, respectively. We then fit the linear low-T 
behavior of these H max (T). Figure 5(b) is a plot of the 
slopes K (black solid circles), obtained from these linear 
fits of the different curves with different percentages of 
MR^™ [Fig. 5(a)], and the corresponding percentage of 
MRf^™ (red solid squares) vs the corresponding values of 
the intercept fields H max (0). We note the sharp increase 
of the slope K, at a certain H max (0) value. 



MATERIALS AND METHODS 

In this Section we provide the supplementary discus- 
sion for the background material that will enable one to 
reproduce our results. First, we outline how the samples 
were prepared and describe the experimental methods 
used to determine the magnetoresistance and its depen- 
dence on magnetic field and temperature. Then, we show 
how describe the procedure we developed to determine 



g-factor and Kondo breakdown 

The origin of the sharp increase in the values of K can 
be interpreted as follows. For a system in the quantum 
critical regime [i.e., low H and T data of Fig. 5(a)], the 
only energy scale is Boltzmann energy £t = k B T. We 
compare this energy scale £t to the quasiparticle Zeeman 
energy, i.e., 

k B T = gp B (H - H QCP ). (1) 
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FIG. 5: (Color online) (a): Temperature dependence of the 
characteristic field H max determined for different percentages 
of the maximum isotropic component of magnetoresistance 
MR^a™ (main panel) and low-T linear behavior of H max for 
< x < 0.20 samples (inset) (b): Slopes K (black solid 
circles) and percentage of MR^" as a function of the corre- 
sponding values of the intercept fields H max {Q). 



SCALING OF THE COHERENCE 
TEMPERATURE 

In this Section, we discuss and offer an interpretation of 
our data for the concentration dependence of the Kondo 
lattice coherence temperature (T^) and superconduct- 
ing critical temperature (T c ). For definiteness, we focus 
primarily on the concentration region x < 0.5, where T cobi 
can be interpreted as a heavy-fermion coherence temper- 
ature of the Ce /-moments. 

We start with the following observation: by combin- 
ing our results on the concentration dependence of the 
coherence temperature and the critical temperature of 
the superconducting transition, we find that to a good 
accuracy both of them are suppressed at the same rate, 
i.e., 

Tll{x) = a Yb T™{x)+(3 Yb , (2) 

where 



From Eq. ([IJ we see that the slope K [Fig. 5(a)] must 
be inverse proportional to the gyromagnetic factor g. So, 
the sharp increase in K is a result of the sharp decrease in 
g. Previous studies 29 have shown that abrupt changes 
in the values of the gyromagnetic factor occur at the 
quantum critical point. Therefore, using this procedure 
we are able to unambiguously determine Hqcp as the 
value of H max (0) at which there is the sharp change in 
the g factor. Specifically, we find that g decreases from 
2.2 just below Hqcp to 1.3 just above Hqcp- 

The experimentally observed changes in the ^-factor 
reflect the transformations that the electronic system 
undergoes under the change in external magnetic field. 
Given that g ~ 2.2 just below Hqcp, numerically close 
to the value found for the free electrons in a metal, sug- 
gests that the conduction electrons become decoupled 
from the local /-moments below the QCP. Thus, this 
observation can be interpreted using the phenomenologi- 
cal theory of "Kondo breakdown" at Hqcp [3D]. Within 
this theory, the changes in the (/-factor are governed by 
the changes in the size of the Fermi surface: larger values 
of the <?-factor correspond to small Fermi surface so that 
the conduction electrons are effectively decoupled from 
the localized /-states. In the opposite limit of smaller 
<?-values, the Fermi surface is large, reflecting the strong 
coupling between the conduction and /-electrons. More 
importantly, the jump in the size of the Fermi surface at 
Hqcp corresponds to the divergence of the quasiparti- 
cle's effective mass. 



a Yb = 16.1 ± 0.2, f3 Yb = 4.09T co/l (x = 0). (3) 

To get better insight into the physical meaning of this 
result, it is instructive to compare Eq. ([2| with the same 
dependence obtained for La-substituted CeCoIn5. We 
find that in this case T co ^ also scales with T c at small 
enough concentration of La [Fig. 2(b) in the main text]: 

T c L a h (x) = a La T c La (x)+(3 La . (4) 

The slight offset between the two curves, Fig. 2(b), ap- 
pears likely because of atomic size differences between 
Yb and La ions. What is truly surprising is that the rate 
it which both T co h and T c are suppressed turns out to be 
the same for both types of substitutions: 

a Yb ^ a La . (5) 

We interpret Eq. ([5| as an indication that the onset 
of the many-body coherence in the Kondo lattice and 
emergence of superconductivity has the same physical 
origin: hybridization between conduction and localized 
/-electron states. In particular, this suggests that Cooper 
pairing develops primarily on the "heavy" (i.e., large) 
Fermi surface. 

Note that, although the suppression rate of Kondo lat- 
tice coherence and superconductivity for both Yb- and 
La-substitutions are the same, both T c ^ and Tj h remain 
much more robust with respect to disorder. This fact can 
be explained by noting that Yb atoms are in a mixed va- 
lence state and, therefore, must be correlated. The cor- 
relation may arise via local lattice deformations. To see 
how the impurity correlations may slow down the sup- 
pression of the coherence temperature, we may consider 
the characteristic length scale R on which impurity dis- 
tribution function significantly deviates from unity. The 
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impurity distribution function determines the probabil- 
ity with which one can find one impurity at a certain 
distance from another. Within the Born approximation, 
one can show that there will be two contributions to the 
self-energy of the conduction electrons. One contribu- 
tion, Ejj, corresponds to the scattering of electrons on 
the same impurity and, upon the averaging over disor- 
der, this contribution is proportional to the concentra- 
tion of impurities n, mp . The second contribution, Sy, 
describes the scattering of electrons on two different im- 
purities and, therefore, is proportional to n? mp . In the 
presence of impurity correlations, however, Sjj becomes 
proportional to nf mp R 3 . Thus, if the radius of correla- 
tions is large enough, ni mp R 3 ~ 1, £.y becomes compa- 
rable with the first, linear in rii mp , self-energy correction 
Consequently, within the large- TV mean field theory, 
one can show that impurity correlations may provide the 
" healing effect" : the rate of suppression of the coherence 
temperature becomes strongly dependent on the impu- 
rity correlation length, R [3"T] . 
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